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Finite- temperature charge-ordering phase transition in quasi one-dimensional (ID) molecular conduc- 
tors is investigated theoretically, based on a quasi ID extended Hubbard model at quarter filling with in- 
terchain Coulomb repulsion V± . The interchain term is treated within mean-field approximation whereas 
the ID fluctuations in the chains are fully taken into account by the bosonization theory. Three regions 
are found depending on how the charge ordered state appears at finite temperature when V± is introduced: 
(i) weak-coupling region where the system transforms from a metal to a charge ordered insulator with 
finite transition temperature at a finite critical value of V±, (ii) an intermediate region where this transition 
occurs by infinitesimal V± due to the stability of inherent ID fluctuation, and (iii) strong-coupling region 
where the charge ordered state is realized already in the purely ID case, of which the transition temper- 
ature becomes finite with infinitesimal V±. Analytical formula for the V± dependence of the transition 
temperature is derived for each region. 
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Charge-ordering (CO) phenomenon has been one of the 
central research subjects in quasi one-dimensional (ID) or- 
ganic 2:1 salts such as (DCNQ^aX'-^ and (TMTTF)2X^*' 
{X: monovalent counter ion). Some members of them having 
a quarter-filled electron or hole band exhibit insulating be- 
havior due to CO. For example, in (DI-DCNQI)2 Ag, an NMR 
measurement clarified the CO phase transition at Tco — 220 
K,' above which non-metallic behavior is akeady seen even 
from room temperature interpreted as due to CO fluctuation.^ 
In (TMTTF)2X, CO is found by NMR as well"* *' and the 
dielectric constant shows a divergence toward Tco. where 
Tco = 70 - 160 K depending on X.^ 

Theoretical studies showed that a minimal model for de- 
scribing such CO is the ID quarter-filled extended Hubbard 
model (EHM) with not only the on-site Coulomb repulsion 
U but also the intersite repulsion ¥7 This model has been 
intensively studied and the ground-state properties have been 
revealed. Using different numerical techniques, the T = 
phase diagram on the iU,V) plane has been obtained with 
high accuracy.** The CO insulator (COI) appears in the large 
{U,V) region. Analytical approaches based on the bosoniza- 
tion theory clarified the mechanism of the transition from 
a Tomonaga-Luttinger liquid (TLL) metallic phase to the 
COI.9-12 

It is useful to consider the TLL parameter Kp depending 
on (U,V) in the metalUc phase of this ID EHM. It takes 
Kp > 1/4 and approaches to Kp = 1/4 at the boundary 
between the TLL and the COI. Based on the TLL theory, CO 
fluctuation is known to develop as xco {T) (x T*^p^^, which 
implies that for 1/2 > Kp > 1/4 it diverges at low temper- 
atures. This divergence does not immediately point to the ap- 
pearance of long range order of CO in the purely ID systems. 
The COI is achieved for larger (U,V), only at T = 0. These 
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features are schematically illustrated in Fig. 1 (a). When U is 
fixed at a large value (e.g., U > At; see Fig. 2 later) and V is 
increased, Kp decreases, and the system transforms through 
three regions bounded hyV = V^^ and V — V^^ defined by 
the values where Kp = 1/2 and 1/4, respectively. The sys- 
tem is the TLL for (i) V < V^^ and (ii) <V< V^'^ 
where in the latter CO fluctuation develops at low-T, and it is 
the COI for (iii) V,}^ < V. 

In order to examine the finite-temperature CO transition 
in the quasi ID materials, in addition to such ID fluctua- 
tion we need to take into account the dimensionality effect by 
the interchain coupling. In these systems two types of inter- 
chain coupling exist in the electronic sector; one is the single- 
particle interchain hopping and the other is the Coulomb in- 
teraction between electrons in the different chains. In this let- 
ter, we consider the latter and discuss the critical tempera- 
ture of CO, Tco. by treating the interchain interaction within 
the mean-field approximation, whereas taking full account of 
the ID fluctuations in the chains by the bosonization theory. 
We will find distinct behavior for the three regions defined 
above : in region (i) the system transforms from TLL metal to 



(a) K^=0 
TLL 



CO^ 




(i) t (ii) t (iii) 

ID y,0 



Fig. 1. Schematic illustration of the finite-temperature phase diagram for 
purely ID case (a) and quasi-lD case with finite V± (b), where the CO 
fluctuation develops in the shaded regions. Here V^^^ express the critical 
value of V for the COI appearing in the quasi-lD case. The quantities V^^ 
and V^^, and the regions, (i), (ii) and (iii) are defined in the text. 
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COI with finite Too at a finite critical value of V±, in (ii) this 
transition occurs by infinitesimal V±, and in (iii) COI in the 
purely ID case always has finite Tco when V± is turned on. 
A schematic phase diagram in the presence of V± is summa- 
rized in Fig. 1 (b). 

OurHamiltonianisgivenbyffqiD = J2j Hlj^+H^.Kere 
Hfj^ represents the j-th isolated extended Hubbard chain: 



and H± expresses the interchain couphng: 



J' 



E 



(1) 



(2) 



to the S = 1/2 XXZ spin-chain coupled by the interchain ex- 
change interaction Jj_(= Vj_) treated in the interchain mean- 
field approach. Several authors have performed such ap- 
proach to the Heisenberg hmit (i.e., V = 2t).'^Bhaseen and 
Tsvehk discussed a nontrivial spectrum of eq. (4) in a param- 
eter region V > 2t,^^ where the system is the COI and has a 
finite excitation gap even for V± = 0.'^ We restrict ourselves 
to a parameter region < V < 2t, where the system is known 
to be metallic if Vj_ = 0, and show that infinitesimaly small 
V± makes the system insulating and gives rise to finite Tco- 
We apply the bosonization method to eq. (4). The den- 
sity operator is given by hi = a{\/2TT)^^dxO{x) — 
7r^^(— l)*sin-\/20(x) with x = ia and a being the lattice 
constant. This parameter region Q < V < 2t corresponds to 
region (ii) with 1/2 > Kp > 1/4 defined above, since it is 
known that in the ID Hubbard model at [/ ^ oo (F = 0), 
Kp takes a universal value Kp 1/2,^° and at ^ = 2t, the 



where ^ is the creation operator of an electron with spin CO transition indicates Kp — 1/4. 



s =T/i at the i-th site in the j-th chain, Uij^s 
l/4andn,.j = n;j.|+nij.|. lneq.(2), denotes the pair 
of adjacent chains and the strength of interchain Coulomb re- 
pulsion is expressed by V± . The number of adjacent chains, z, 
for (DCNQ1)2X and (TMTTF)2X would be z = 4 and z = 2 
judging from the crystal structures, although their interchain 
networks are rather complicated, as we will discuss later. We 
do not take into account the interchain hopping which should 
not yield quahtative changes to the present analysis at the 
temperature above the crossover energy scale due to the in- 
terchain hopping.'^ 

We treat eq. (2) in the interchain mean-field approach, 
which is known to be effective in the weak V± region, by 
considering the Wigner-crystal-type CO pattern stabilized in 
the ID EHM. We assume the CO pattern to be anti-phase be- 
tween adjacent chains, which is naturally expected to gain the 
Vj_ term. The resulting effective ID Hamiltonian is written as 

+ ^ E '^''T '^»4 + ^ E '^'+1 



(3) 



where the chain index j is suppressed and N is the total num- 
ber of sites in a chain. The amphtude of CO is written as n. 

To obtain a qualitative insight of this model we first discuss 
the U oc limit, where some exact results for the ID model 
based on the Bethe ansatz can be used. In this case, we can 
neglect the spin degree of freedom which acts freely, and the 
charge degree of freedom is reduced to a half-filled spinless 
fermion model: 

Hu^oo = -t ^((ijc^i+i + li-c-) +V'^hi hi+i 

i i 

+ zV^nY.i-'^yni + zNV±ny2, (4) 

i 

where the creation operator of the spinless fermion is ex- 
pressed by dl and hi = d\di — 1/2. This model is equivalent 



In the bosonized form of the Hamiltonian, there appears 
a nonhnear term, y/(27r^o) cos2v^^, which becomes ir- 
relevant but has an effect to renormalize the velocity, the 
TLL parameter and the prefactor of the bosonized density 
operator. By taking into account these effects properly, the 
density operator is given by h, = a(V2^)-ia,0(.x) - 
C7r~^(— 1)^ sin-\/2^(a;), where c is a nonuniversal constant 
depending on V/ 1 and its exact value has been suggested.^' 
We will use the renormalized velocity and the TLL parame- 
ter, and neglect the cos 2-J29 term. After this procedure, we 
switch on the V±_ term, and then the effective Hamiltonian 
density reads 



V 



5-L 



■{d,0f + 2Kp{d^ 



sm\/2e + —Vi^v? , 
2a 



(5) 



where v = nta^Jl - {V/2tf / aiccos{V/2t), Kp = [A- 
(4/7r) cos-i(y/2t)]-\ and = 27TCzV±an. The TLL pa- 
rameter Kp determines the nonuniversal constant c = c(Kp). 

Since the g± term in eq. (5) is a relevant perturbation, when 
n 7^ the system is not the TLL but has an excitation gap,^^ 

^ l/(2-2A:p) 

9± 



m{n) 



(6) 



where = 7r-i(7r/4)i/(P+i)r(p/(p+l)) T-^{l/{p+l)) 



[r((p-hl)/2) T-^{p/2)] 
per site is given by^^ 



2/(p+l) 



. The ground- state energy 



E(n) = -—m^(n) tan [ - ^ 



2K, 



zV±n'^/2. (7) 



The CO amplitude n is determined from the condition to min- 
imize eq. (7). The optimized magnitude n is given by 



tan(7rXp/(2 - 2Kp)) 



2-2Ko 



<Kp) 



2-kk{j 



T/ \ Kp/(l-2Kp) 



2v 



(8) 
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Inserting this expression to eqs. (6) and (7), the excitation 
gap and the ground-state energy are given by m = m(n) oc 

yl/(2-iK,) ^ ^ ^^^^ ^ yl/{l-2K,)_ ^^^^^ 

finitesimal value of interchain interaction transforms the sys- 
tem from TLL in the ID case to the COI for all values of 
< V < 2t. Tco can also be estimated in a similar way. The 
free energy up to is given by 



fin) 



6v 



+ 



2a 



(9) 



where b{Kp) 



'{Kp) sin{2TT Kp)B^{Kp,l - 2Kp) with 



B{x,y) = r{x)r{y)/r{x + y). From the condition where 
the coefficient of in eq. (9) vanishes, Tco is determined as 

'-V^a 



Tco 



2'Ka 



zKKp), 



2'kH 



(10) 



KxV = 2t, Tco is proportional to Vj_ . " In this case, the ra- 
tio of the gap at T = and the transition temperature becomes 
'm/Tco\v=2t- ~ 2.47, which reproduces the result obtained 
in Ref. 23. The formulas of the amplitude n [eq. (8)] and the 
transition temperature Tco [eq- (10)] are valid for V<V^^, 
since these diverge when V ^ due to the unphysical ultra- 
violet divergence in eqs. (7) and (9). This problem could be 
resolved simply by introducing an ultraviolet cutoff. 

Now we consider the case of finite U. An effective Hamil- 
tonian for low-energy states (k « ±7r/(4a)) can be obtained 
by integrating out high-energy states.'^ The Hamiltonian is 
separated into the charge part describing the CO transition and 
the spin part. The latter is essentially the same as the Hamil- 
tonian of ID isotropic Heisenberg model and leads to gapless 
excitations both in the TLL and the COI. The density operator 
can be bosonized as 

rii =—dxOp{x) + cos{2kYX + 9p{x))sm6„{x) 
-|-C2^^^cos2^p(a;), 



(11) 



where 9p{x) and 6^ (x) are the bosonic phase variables for the 
charge and spin degrees of freedom. c\ and C2 are nonuniver- 
sal numerical constants. In the absence of n, the charge part 
is written by the following phase Hamiltonian, 

1 



47r 



-{djpf ^KpQ{dx<t>pf 



cos46>p, 

2(7roj^ 

(12) 



where [9p{x),—dx'4ip{x')/{2TT)\ — i5{x — x'). The charge 
velocity Vp, the parameter Kpo, and the magnitude of the 1/4- 
filled Umklapp scattering are nonuniversal and its ex- 
pressions obtained perturbatively are given in Ref. 12. In the 
presence of the interchain mean-field n, the additional Umk- 
lapp scattering appears, as in the last two terms in eq. (5), 



n' = —V±ncos20p + ^V±n^. 



(13) 



By minimizing {Hp + H') with respect to n, the amplitude n 
is determined as 



— (cos2(?,, 



(14) 



where (• • • ) denotes the finite-T expectation value with re- 
spect to Hp + H'. We note that the parameter C2 is propor- 
tional to the on-site interaction, C2 ~ U/{-j2iTt) within the 



perturbative treatment. This implies that the CO fluctuation 
with 2-fold periodicity can emerge due to the correlation ef- 
fect. 

In order to estimate the second-order transition tempera- 
ture, the lowest-order perturbation theory in n is sufficient. 
By expanding the r.h.s. of eq. (14) with respect to H' , the 
equation to detemune Tco is obtained as^^ 

r-l/Tco 



1 



dx 



dr 



X {TrCos2ep{x,T)cos20p{O,O))Q\T=Tao, (15) 



where (• • • )o is the expectation value with respect to 
Hp. In the following, we employ naive renormalization 
group (RG) arguments to estimate Tco- By neglecting 

the anisotropy in space and time, the correlation function 
(Tr cos26p{x, t) cos 29p{0, 0))o can be roughly calculated as 

{Tr COS 29p{x, t) cos 29p{0, 0))o 



■ exp 



ln[min(rp,{p)/a] 



dl{4Kp{l) + 2Gi/i{l)) 



(16) 



where = sjx^ + {^p^Y ^i^^ ~ '"pH'^T) is the char- 
acteristic thermal coherence length. The quantities Kp(l) and 
Gi/i{l) are the solution of the RG equations, 

dKp{l)/dl = -8Gl,S)Kl{l), (17a) 

dGi/4(/)M = (2 - d,Kp{l))G^,S), (17b) 

where the initial conditions are Kp{Q) = Kp^ and 
^1/4(0) = fl'i/4/(27rt'p). These RG equations have 
two kinds of fixed points for realistic parameters. One 
is {Kp{oo),Gi/4{oo)) = (0, — 00) and the other is 
{Kp{oo), 01/4(00)) = {Kp, 0). For VI = the former cor- 
responds to the COI at T = 0, while the latter to the metallic 
TLL state. From eqs. (15) and (16), Tco can be estimated by 
F (bvp/Tco) = 1 where 



_zc\aV±_ 



TTVo 



In s 



X exp 



dy 



dl{4Kp{l) + 2Gi/4{l)-2) 



(18) 



and 6 is a positive numerical constant. 

Concerning the realization of the COI at finite-T, we find 
three distinct regions on the iU,V) plane shown in Fig.2. In 
region (iii) where the COI is already realized in the purely ID 
case, infinitesimal V± makes Tco finite, as is naturally ex- 
pected. In addition, there is a region (ii) with 1/2 > Kp > 



10 



U/t 




(iii) 



(i) ^ *:p=o.5 



10 



V/t 



Fig. 2. Three distinct regions classified by appearance of the COI at the 
finite temperature. 
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Table I. The V± dependence of the transition temperature Tqq in the re- 
spective regions, where A, B and C are positive constants. Corresponding 
Kp values are also listed. In the region (i), the transition happens only for 
V± > A{4Kp - 2). 







Tco 


(i) 
(ii) 
(iii) 


Kp > 1/2 

1/4 < Kp< 1/2 

not defined 


(X [1 - A{iKp - 2)/yL]i/(4Kp-^; 
oc [1 - A{4Kp - 2)/V±]^/('^'^p-'^"l 
= V±[l + Bln{C/V±)] 



1/4 where infinitesimal V± also gives rise to finite Tco even 
in the metalUc region for V± = 0, since F{oo) diverges due 
to the divergence of CO fluctuation while Gi/4,{oo) = 0. This 
is consistent with the case of [/ ^ oo discussed above, where 
< V < 2t correspond to region (ii) in this limit, as men- 
tioned. Moreover, even in the region (i) without divergence of 
CO fluctuation, finite amount of V± makes the system COL 
The obtained expressions for Tco as a function of V± is sum- 
marized in Table. I. The critical value of V± for appearance of 
the COI in the region (i) is proportional to AKp — 2. We note 
fliat Tco oc exp{-A/V±) at Kp = 1/2, and Tco oc V± for 
Kp = 1/4, also consistent with the f7 — > oo case. 

Now let us discuss the relevance of our results to the exper- 
iments. In (DI-DCNQI)2Ag, the non-metallic behavior above 
Tco - 220 might be due to the 1/4-fiUed Umklapp scat- 
tering described in eq. (12). The peak at T = Tco in the 
derivative of the resistivity as a function of T in Ref. 3 seems 
to be triggered by the generation of additional Umklapp scat- 
tering cos 29 p in eq. (13) due to the appearance of the CO. In 
(TMTTF)2X, the system is also insulating already above Tco 
and a change in the slope of flie resistivity curve is observed.-^* 
However, in this case the slight dimerization in the chain di- 
rection gives rise to another relevant 1 /2-filled Umklapp scat- 
tering.^^ Whether or not tiie 1/4- filled Umklapp scattering is 
effective above Tco is difficult to judge, but the change in the 
slope should be due to the additional Umklapp scattering by 
CO as discussed in this paper. 

In these quasi ID compounds, the interchain Coulomb in- 
teraction is in fact considered to be fairly large. For exam- 
ple in some members of (TMTTF)2X it is even estimated to 
be comparable to the intrachain ones in quantum chemistry 
calculations.^^ This is noticeable since the interchain elec- 
tron hopping is one order of magnitude smaller than that in 
the intrachain direction due to the anisotropic shape of the 
HOMO of the TMTTF molecule. However we should note 
also that the actual interchain networks are rather compli- 
cated. In the structure of (DCNQI)2X, it has a spiral symme- 
try when rounding each plaquette in the network. It is pointed 
out that this gives rise to frustration for the Wigner-crystal- 
type CO since the periodicity of the spiral does not fit the 2- 
fold period along the DCNQl chains.^*- In (TMTTF)2X, such 
frustration is more straightforwardly expected since the in- 
terchain coupling is zigzag-like therefore anti-phase and in- 
phase CO patterns would have close energy. These should re- 
duce the effective interchain interaction to a smaller value at 
a first approximation. 

In the present analysis, we considered the Wigner-crystal- 
type CO with a 2-fold periodicity only, namely, the "4fcF" 
CDW state. In the small V region, however, the 2fcF CDW 
state having a 4-fold periodicity may also be stabilized, since 
its fluctuation develops as xf"^ {T) oc T^"'^. From the sim- 
ple power counting, one finds that 2kF fluctuation becomes 
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dominant if Kp > 1/3. Thus in the phase diagram where 
1/3 < Kp < 1/2, competition and/or coexistence of the CO 
and the 2kF CDW may be possible. This problem needs fur- 
ther investigation. 

In conclusion, we investigated the finite-temperature CO 
transition in the quasi- ID electron system at quarter filling 
coupled by the interchain Coulomb repulsion V±, by using the 
interchain mean-field theory and the bosonization method. It 
was shown that the interchain interaction gives rise to finite 
transition temperature Tco, with different behavior depend- 
ing on the parameter in the ID limit. When it is in the COI, 
V± gives rise to finite temperature phase transition. In the TLL 
phase, in the parameter region where the charge-fluctuation 
develops, i.e., 1/2 > Kp > 1/4, infinitesimally small V± 
also produces the COI with finite Tco, while for Kp < 1/2, 
at a finite critical value of V± the system transform from the 
TLL to the COI. 
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